University  of  Missouri-Columbia 


Weibull  Accelerated  Life  Tests  When 
There  are  Competing  Causes  of  Failure 


John  P.  Klein 

Department  of  Statistics,  Ohio  State  University 
Asit  R  Basu 

Department  of  Statistics,  University  of  Missouri-Columbia 


Technical  Report  104 
Department  of  Statistics 


June  1961 


Mathematical 


Sciences 


<0^ 


81  8  06  039 


SECURITY  CLASSIFICATION  of  this  RACE  (Whan  Data  F.ntarad ) 


REPORT  DOCUMENTATION  PAGE 


I.  REPORT  NUMBER 


«  T|TLE1a««  Sub, Ilia) 


a  Arp  HEAD  INSTRUCTIONS 

_ -BpFORE  COMPLETING  FORM 

|2.  GOVT  ACCESSION  NO.  l^IEClPfENT'S  CATALOG  NUMBER 


Weibull  Accelerated  Life  Tests  When  There  Are 
Competing  Causes  of  Failure  _ . . . . 


Technical 


Report 


4.  PERFORMING  ORG.  REPORT  NUMBER 


N0P014-78  C  J9655 


auxuoaoi  ^ — j  .».lqmtalact  o*-vr*nt  voCiERfij 

John  P./ Klein  (/£  N0P014-78  C  J9655 

As  it  P./Basu  /  ^  ‘  ^ 


PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Department  of  Statistics 
University  of  Missouri-Columbia 
Columbia,  MO.  65211 


it.  controlling  office  name  ano  address 
Office  of  Naval  Research 
Department  of  the  Navy 
Arlington,  VA.  22217 


14.  MONITORING  AGENCY  NAME  «  AOOREESrtl  dlllatanl  ltam-C*Uiaiflnl  Olll ca)  IS.  SECURITY  CLASS,  (ol  thla  raport) 


ADORES $(l^dyler. 


IJU 


Unclassified 


IS*.  DECL  ASSIFIC  ATION/  DOWNGRADING 
SCHEDULE 


ENT  (ol  the  abetted  entered  In  Block  20 ,  II  dllloront  Irom  Rmport) 


it.  KEY  WORDS  (Continue  on  reveree  aide  II  neceeeary  and  Identity  by  block  number) 

safe  dose  levels,  competing  risks,  accelerated  lifetests, 
Hartley  and  Sielkin  model. 


IB.  ABSTRACT  (Continue  an  reeeree  aide  II  neceeeaey  and  Identity  by  block  number) 

Accelerated  life  testing  of  a  product  under  more  severe  than  normal  condi¬ 
tions  is  commonly  used  to  reduce  test  time  and  costs.  Data  collected  at  such 
accelerated  conditions  are  used  to  obtain  estimates  of  the  parameters  of  a 
stress  translation  function.  This  function  is  then  used  to  make  inference 
about  the  product's  life  wider  normal  operating  conditions. 

We  consider  the  problem  of  accelerated  life  tests  when  the  product  of  inter¬ 
est  is  a  p  component  series  system.  Each  of  the  components  is  assumed  to  have 
an  independent  Weibull  time  to  failure  distribution  with  different  shape  para- 


20 i  -^meters  and  different  scale  parameters  which  are  increasing  functions 
of  the  stress.  A  general  model  is  used  for  the  scale  parameter  which 
includes  the  standard  engineering  models  as  special  cases.  This 
model  also  has  an  appealing  biological  interpretation. 

Maximum  likelihood  estimators  of  the  component  parameters  are  ob¬ 
tained  for  type  I,  type  II,  and  progressive  censoring.  These  esti¬ 
mators  are  used  to  obtain  estimators  of  the  probability  of  component 
and  system  survival  under  normal  operating  conditions.  This  method 
is  illustrated  by  an  example. 


/ 


in  lifcH 


•  •  if  •■Hpw  4rm 


WE I BULL  ACCELERATED  LIFE  TESTS  WHEN  THERE  ARE 
COMPETING  CAUSES  OF  FAILURE 


John  P.  Klein 


Department  of  Statistics 
The  Ohio  State  University 
Columbus,  Ohio  43210 


Asit  P.  Basu 

Department  of  Statistics 
University  of  Missouri 
Columbia,  Missouri  65211 


Key  Words  and  Phrases:  safe  dose  levels,  competing  risks, 
accelerated  lifetests,  Hartley  and  Sielkin  model. 


ABSTRACT 

Accelerated  life  testing  of  a  product  under  more  severe 
than  normal  conditions  is  commonly  used  to  reduce  test  time 
and  costs.  Data  collected  at  such  accelerated  conditions  are 
used  to  obtain  estimates  of  the  parameters  of  a  stress  transla¬ 
tion  function.  This  function  is  then  used  to  make  inference 
about  the  product's  life  under  normal  operating  conditions. 

We  consider  the  problem  of  accelerated  life  tests  when 
the  product  of  Interest  is  a  p  component  series  system.  Each 
of  the  components  Is  assumed  to  have  an  Independent  Wei  bull 
time  to  failure  distribution  with  different  shape  parameters 
and  different  scale  parameters  which  are  increasing  functions 


of  the  stress.  A  general  model  is  used  for  the  scale  parameter 
which  includes  the  standard  engineering  models  as  special 
cases.  This  model  also  has  an  appealing  biological  interpreta¬ 
tion. 

Maximum  likelihood  estimators  of  the  component  parameters 
are  obtained  for  type  I,  type  II,  and  progressive  censoring. 
These  estimators  are  used  to  obtain  estimators  of  the  probab¬ 
ility  of  component  and  system  survival  under  normal  operating 
conditions.  This  method  is  illustrated  by  an  example. 


1.  INTRODUCTION 

Accelerated  life  testing  of  a  product  is  often  used  to 
obtain  information  on  its  performance  under  normal  use  condi¬ 
tions.  Such  testing  involves  subjecting  test  items  to  condi¬ 
tions  more  severe  than  encountered  in  the  item's  everyday  use. 
This  results  in  decreasing  the  item's  mean  life  and  leads  to 
shorter  test  times  and  reduced  experimental  costs.  In  engi¬ 
neering  applications  accelerated  conditions  are  produced  by 
testing  items  at  higher  than  normal  temperature,  voltage, 
pressure,  load,  etc.  In  biological  applications  accelerated 
conditions  arise  when  large  doses  of  a  chemical  or  radiological 
agent  are  given.  In  both  cases  the  data  collected  at  the  high 
stresses  is  used  to  extrapolate  to  some  low  stress  level  where 
testing  is  not  feasible. 

Several  authors  have  considered  the  problem  of  analyzing 
accelerated  life  tests  when  the  product  has  only  a  single  mode 
of  failure.  Nelson  (1974a)  has  a  bibliography  of  applications 
and  analysis  of  each  tests.  Mann,  Schafer,  and  Singpurwalla 
(1974)  derive  least  squares  and  maximum  likelihood  estimators 
of  model  parameters  when  the  underlying  failure  distribution  is 
exponential.  Nelson  (1972c)  describes  a  graphical  solution  to 
this  problem. 


Nelson  (1970)  derives  graphical,  maximum  likelihood  and 
least  squares  estimators  of  model  parameters  when  the  under¬ 
lying  distribution  is  Wei  bull.  Meeker  and  Nelson  (1974a  and 
1974b)  derive  maximum  likelihood  estimators  of  the  model  para¬ 
meters  in  the  Wei  bull  case  when  the  data  is  type  I  or  type  II 
censored.  They  also  discuss  the  optimal  strategy  for  designing 
such  tests.  Mann  (1972)  has  also  discussed  optimal  design 
strategy.  Tolerance  bounds  for  the  Weibull  model  are  discussed 
in  Mann  (1978). 

Several  papers  have  been  written  on  analyzing  accelerated 
life  tests  when  the  failure  distribution  is  normal  or  log 
normal.  In  a  series  of  papers  Nelson  (1971,  1972a  and  1972b) 
has  considered  maximum  likelihood,  least  squares,  and  graphical 
estimation  procedures  for  an  Arrhenius  model  when  all  failure 
times  are  known.  For  this  model  it  is  assumed  that  mean  of  the 
log  failure  time  is  linear  in  the  stress,  and  that  the  variance 
is  independent  of  the  stress.  Nelson  and  Hahn  (1972,  1973) 
derive  best  linear  unbiased  estimators  of  the  regression  para¬ 
meters  of  this  model  for  type  II  censored  samples.  Kielpinski 
and  Nelson  (1975)  discuss  maximum  likelihood  estimation  pro¬ 
cedures  for  this  model  when  the  sample  is  type  I  censored. 

Several  papers  have  been  written  on  analyzing  accelerated 
life  tests  when  more  than  one  failure  mode  is  present.  Here 
failures  can  occur  from  any  one  of  k  independent  causes.  Sam¬ 
ple  information  consists  of  a  failure  time  for  each  item  and 
the  cause  of  failure.  Assuming  that  for  a  given  stress  V  each 
failure  mode  follows  an  independent  log  normal  distribution 
with  parameters  u^(V)  =  V,  and  o^2  constant  with 

respect  to  V,  i  =  1,  ...,  k,  Nelson  (1973)  obtains  graphical 
estimates  of  and  6^  when  there  is  no  censoring.  For 

this  model,  Nelson  (1974b)  obtains  maximum  likelihood  esti¬ 
mates  of  and  2 . 


Klein  and  Basu  (1980a)  have  considered  the  above  problem 
when  the  component  lifetimes  are  exponentially  distributed  and 
the  data  is  type  I,  type  II  or  progressively  censored.  Klein 
and  Basu  (1980b)  have  described  the  analysis  of  accelerated  life 
tests  of  series  systems  when  the  component  failure  times  follow 
a  Weibull  distribution  with  common  shape  parameter.  The  object 
of  this  paper  is  to  consider  the  case  when  the  shape  parameters 
are  different. 

In  section  2  we  present  the  model  to  be  used  for  acceler¬ 
ated  life  tests  in  the  competing  risks  framework.  In  section 
3  we  use  this  model  to  analyze  accelerated  life  tests  where 
there  are  competing  causes  of  failure  and  the  data  is  type  I, 
type  II,  or  progressively  censored.  Finally,  in  section  4  an 
example  is  presented. 

2.  THE  MODEL 

The  problem  considered  in  the  sequel  is  as  follows.  Consi¬ 
der  a  p  component  system  with  component  lifetimes  X1 ,  X2,  ....  Xp. 
Suppose  that  under  normal  stress  conditions  these  components  have 
long  lifetimes  making  testing  at  such  conditions  unfeasible.  To 
reduce  test  time  and  cost,  s  stresses,  ,  . ..,  V$  are  selected 
and  a  life  test  is  conducted  at  constant  application  of  the 
selected  stress.  We  wish  to  use  this  information  to  make  infer¬ 
ence  about  the  component  lifetimes  under  normal  stress  conditions. 

Consider  the  following  model  introduced  by  Klein  and  Basu 
(1980b)  elsewhere. 

At  a  stress  V^,  i  =  1,  ....  s  assume  that  the  jth  component 
has  a  hazard  rate  given  by 

hj(x.  v,;  SLj.  Ij)  *  9j(x.  SjIVV  (2J) 

"I  s  «  i  S  J  *  1$  *  •  •  i  p« 

Oj)  a  Weibull  form  is  assumed,  that  is 


9J<X-  °J>  =  “j1 


(2.2) 


,  Oj  >_  0,  t  >  0. 

The  ou ' s  may  vary  from  component  to  component  to  allow  for  dif¬ 
ferences  in  component  reliability. 

For  A.j(V,  B-)  we  assume  a  model  of  the  form 

J  J 

kj 

\j(v.  Ij)  =  e*p  (  E  Bj£ejA(v) )  •  (2-3) 


where  0jO(v)  =  1  and  8j-|(V),  •••»  ^.(V)  are  kj  nor,-decreasing 

J 

functions  of  V.  The  0 - ( * ) ' s  may  differ  from  one  component  to 

j 

another. 

This  model  includes  the  standard  models,  namely,  the  power 

6., 

rule  with  Aj(V,  J  ;  the  Arrhenius  reaction  rate  model 

with  Aj(V,  Bj)  =  exp(Bj0  -  Bj^/V);  and  the  Eyring  model  for  a 

Bi-i 

single  stress  with  Aj(V,  Bj)  =  V  exp(Bj0  -  Bj2/V)  as  special 


cases. 

The  model  also  can  be  derived  from  the  interpretation  of 
the  effects  of  a  carcinogen  on  a  cell  as  proposed  by  Armitage 
and  Doll  (1961).  For  details  see  Klein  and  Basu  (1980b).  To 
produce  cancer  in  a  single  cell,  k  independent  events  must 
occur.  The  effect  of  an  increased  dose  of  a  carcinogen  is  to 
increase  the  rate  at  which  these  k  events  occur.  If,  for  the 
j**1  disease,  this  increase  is  of  the  form  exP(8j£9j£(v))  for 
£  ■  1,  ...»  k.  the  model  (2.3)  is  obtained.  If  this  increase 
is  assumed  linear  the  model  of  Hartley  and  Sielkin  (1977)  is 
obtained.  Thus  the  model  of  Hartley  and  Sielkin  is  a  first 
order  Taylor  Series  approximation  to  (2.3)  when  0j ^ ( V)  =  V  for 
l  1  li  •••»  ^ j « 

Consider  an  accelerated  life  test  conducted  at  constant 
applications  of  s  stress  level,  V^,  ...,  V$.  Let  Xil  ’  Xi2* 
...,  denote  the  component  lifetimes  of  the  p  component 


series  system  put  on  test  at  stress  V^.  Assume  that  the  compo¬ 
nent  lifetimes  are  independent.  We  are  not  allowed  to  observe 
X 

minimum(X^i 

which  of  the  p  components  is  the  minimum. 


,  ....  Xip  directly  but,  instead,  we  observe  Y^  = 

X^)  and  an  indicator  variable  which  describes 

We  shall  use  the 


method  of  maximum  likelihood  to  estimate  cij  and 


(B 


JO' 


)»  j  =  1,  ....  P  for  various  censoring  schemes. 
J3.  ESTIMATION  OF  PARAMETERS 


3.1  Type  I  censoring 

For  this  censoring  scheme  n..  items  are  put  on  test  at 

stress  V^,  i  =  1,  ...»  s.  The  £*h  system  on  test  at  stress 

is  tested  until  it  fails  or  until  some  fixed  time  x^  at  which 

it  is  removed  from  the  study.  The  x^'s  may  very  from  item  to 

item  to  allow  for  staggered  entry  into  the  study.  Let  r^  be  the 

number  of  systems  which  fail  prior  to  their  censoring  time  at 

stress  V^.  Let  r^  denote  the  number  of  these  whose  failure  was 
1  1  ^  th 

caused  by  failure  of  the  j  component.  Let  X^  denote  the 
failure  time  of  the  r-j  systems  whose  failure  is  due  to  failure 
of  component  j.  For  convenience  let  Y^,  i  =  1,  ...»  s, 
i  a  1,  ....  r..  denote  the  failure  times  of  the  r^  systems 
regardless  of  the  cause  of  failure. 

The  overall  log  likelihood  can  be  written  as 


P 

inL  =  l  tnL. 
j=l  3 


(3.1.1) 


where 


"  wexpy0  j  =  1  •  •••• 


p 


(3.1.2) 


where 


1 


s  (3.1.3) 


with  i  =  1,  . ..,  s,  a  =  1,  ....  n^-r^  the  censoring  times 
of  the  n^-r^  systems  removed  from  the  accelerated  life  test. 
When  all  items  on  test  at  stress  have  a  common  censoring 


1  i 

time,  Tj,  then  T^otj)  =  l  Yi  J 

A”  1 


The  likelihood  equations,  which  must  be  solved  numerically 

^  ✓N  A  A 


for  the  maximum  likelihood  estimators, 
are: 


aj*  Bj0*  Bjl» 


6  Jin 

0  =  - j- = 

oa . 

J 


kj 

'  Till>(“j)«P  (J0  BjtVV  ) 


j  ■  1. 


(3.1.4) 


where 


Ti  =  Yifc*nYU  +  TUdtnTiV  1  =  1»  s, 

(3.1.5) 

and 

6£nL  s  kj 

°  ill  rUeJ«(V  -  Ti(“j>eju<vi>«P  (J0  eJteJs.(Vi))  . 

j  =  1  *  •  ••»  P  U  ~  0,  ...,  kj  (3.1.6) 


| 

i 

j 


The  second  partial  derivatives  of  the  log  likelihood  are 


p.  (3.1.7) 


and 


ni_ri 

^•j  Tii^nTi£,)2»  =  •••*  P» 


<52*nL.  s  .  . 

=  j,  xijejU(''i>Ti0)(“j) 

J  —  1 1  ■  ■  • »  p  u  =  0|  . . .  i  k . 

J 


'it 


(3.1,8) 


(3.1.9) 


62inLj  s 

«0juejw  =  1|1  XijTi(aj)VVi)6ju(Vi) 

w  1*  ...j  k .  u  =  Of  . , .  j  k*. 

J  j 

To  find  the  information  matrix  let 

1  ,f  YU  iT1f  1  =  1 . s.  1  *  >.  • 

0  otherwise 

and  define 


(3.1.10) 


» 


(3.1.11) 


*  p(CU  =  1)  =  1  -  exp  ( 


1  s  1 . s  l  =  1 . n^  (3.1.12) 

The  conditional  density  function  of  Y.  given  C.ff  =  1  is 

1  ~  t  £ 


'kill'll  ■  ’>  ■ 


P  a,-l 

l  “jx1j* 


ill 


a . 


‘iS, 


‘  0  otherwise. 


exp  (-  X  Xijy  J)  .  y  <  TU 


i=l 


(3.1.13) 


Now 


E(T,(aj»  -  E  ( j|  -Cu)t“{) 


-1  l  ^  (  Jl  %Ximy  6XP  m-1  Ximy°"  ]  ^ 

+  1  TU  exp  (  -  1  XimTu) 


i  1 )  *  •  •  j  P  3  "  1 »  • • •  i  s. 
where  the  integral  must  be  evaluated  numerically.  Similarly, 

E(Ti(1)(°j»  = 


(3.1.14) 


l  l  ny  (  f  Vim* m  )  exP^  -  1  ximya'n) dy* 

£=1  0  m=1  m  inr  m=l  inr 

+  T1J(Anxu)exp(  X^J).  (3.1.15) 


j 

i 

i 

i 

i 

t 


and 


E(T(Z)(aj))  • 


l  f'U y°W  (  l  Vl^V’ )  «P  (  -  t  **/'' 

S,=l  o  m=l  m=l 


ni  a. 


+  J,  Tii(fi'nTi£)2exp(  J,  AijTU  ) 


i  =  1 . s  j  =  1 , 


(3.1.16) 


Now  the  probability  that  the  system  fails  prior  to  time 
due  to  failure  of  cause  j  at  stress  i  for  i  s  1,  ...»  s, 
j  *  1,  ...,  Pi  i  s  1,  ...»  n.  is 


;T%uVe*p(-  l 

0  J  m=l 


X.  u  m)du 


(3.1.17) 


which  must  be  evaluated  numerically.  Hence 

ni 

E(rjf)  =  l 


(3.1.18) 


The  asymptotic  covariance  matrix  can  now  be  obtained  by 
using  (3.1.18),  (3.1.16),  (3.1.15),  and  (3.1.14)  to  calculate 
the  expected  values  of  (3.1.7),  (3.1.9),  and  (3.1.10).  An 
estimator  of  this  matrix  can  be  obtained  by  substituting  a-,  and 

a  •  ^ 

In  the  appropriate  expressions. 

J 

3.2  Type  II  censoring. 

For  this  censoring  scheme  n^  systems  are  put  on  test  at  each 
of  the  s  stress  levels  and  testing  continues  until  a  preassigned 
number  r^  have  failed  at  which  time  testing  is  stopped.  Suppose 
that  r^j  systems  fail  due  to  failure  of  the  jth  component, 
j  «  1,  ...,  p.  Let  ,  ....  X-jjr^  denote  the  failure  time 
of  those  r,-j  systems  at  stress  Vi  whose  failure  was  caused  by 


'’WWPT' 


r 


failure  of  the  jth  component,  1  =  1,  ....  s,  j  =  1,  ....  p, 
l  -  1,  ....  r^..  Let  Y^ij,  ....  denote  the  ordered 

failure  times  of  the  r^  systems  observed  to  fail  at  stress  1, 
regardless  of  the  mode  of  failure. 

One  can  show  that  the  likelihood  of  interest  is  given  by 
(3.1.1)  and  (3.1,2)  wtth. 


i  *  1 . s  j=l»...»p.  (3.2.1) 

The  likelihood  equations  are  given  by  (3,1.4)  and  (3,1,5)  with 


T1°S>  ■  j’  +  ("1  -  '-i>Yiir1)ilnVi(r1) 


i  8  1,  s  j  =  1 ,  . . . ,  p.  (3.2.2) 

The  matrix  of  second  partial  derivatives  is  given  by  (3.1.7)  to 
(3.1.9)  with 


v> 


(“j>  •  J,  Yl|t)(^1U)l2  *  ("1  - 

j  =  1 . p  1  8  1,  ...»  s.  (3.2.3) 

To  find  the  information  matrix  note  that  the  density  of 
is 


a,-l 


- — ^ -  (?  }  1  -  exp  (-  f  X^yJ.x ) 

(H  -  l).1  (ni  -  i)\  b-1  J  lW  1J  1111 

-(ni  -1  +  1)(k4)) 

J  * 


i-1 


l  »  1 


9  •  •  •  9 


n.j  1  =  1 ,  . . . , 


s. 


(3.2.4) 


Ed^oj))  = 


l  «>  a. 


J,  l  + 


-  r()  ||  ylfrj)f(y1(r1)'dyt(r1)' 

E<T11)(aj>)  =  J,  /Q  ylll)tny1(t)f<yiU)>dylU) 
*  (n1  -  ri>  / 

E(T!2)(“j»  -  J,  Jy3t)<*"y1U)>’f<yiU)>dyiU) 


+  (ni  -  r^ 


^  —  •  •  •  i  S  3  ■  1  j  •  *  ■  i  p  • 


These  integrals  must  be  evaluated  numerically.  Now  with  r. 
fixed,  (r^,  ....  r^)  has  a  multinomial  distribution  with 
parameters  given 

"ij  *  /"  xu“j“aj',exP  ( -  Jm/1)  du- 

n  • 


E(rij)  =  r1nij* 

The  information  matrix  can  be  obtained  by  using  (3.2.5), 
(3.2.6),  (3.2.7),  and  (3.2.8)  to  calculate  the  expected  values 
of  (3.1.7),  (3.1.9),  and  (3.2.10).  An  estimate  of  this  matrix 

A  A 

can  be  obtained  by  substituting  and  8j  in  the  appropriate 
expressions. 


1 


3.3  Progressive  censoring 

For  this  censoring  scheme  N,  items  are  put  on  test  at  the 

+  h  * 

itn  stress  level.  Let  -rn,  ...»  be  fixed  censoring  times 


at  which  a  fixed  numer  of  items*  c^,  ....  c^  are  removed  from 

the  test.  At  time  either  a  fixed  number  c^  items  are 

removed  from  the  test  or  testing  is  terminated  with  a  random 
number,  c^  ,  systems  still  functioning.  Assume  that  is 

sufficiently  large  to  allow  removal  of  the  required  number  of 
items. 

Mi 

Let  nf  s  Nj  ■  y  c.-K  be  the  number  of  items  which  are 
i  i  ik 

observed  to  fail  and  let  Y^,  it  =  1 . n^  denote  their  failure 

times.  That  is,  Y^,  ...,  Y^n  are  the  n^  system  failure  times 

regardless  of  the  mode  of  failure.  Suppose  that  r. .  of  the  n^ 

t’h  ^  ' 

failures  were  caused  by  failure  of  the  j  component,  and  the 

respective  failure  times  are  X^,  ....  X1-jr  ,  i  =  1 . s, 

j  *  1,  ...»  p. 

The  likelihood  function  of  interest  is  given  by  (3.1.1) 
which  can  be  factored  into  component  likelihoods  as  in  (3.1.2) 
with 

ni  a.  ^i  a. 

Ti^aj^  =  Yiit  +  ^  TiJtciJt’  j  *  1 . p,  i  =  1,  ...»  s. 

The  likelihood  equations  are  as  in  (3.1.4),  (3.1.5)  with 


rij  ^ 

=  JiYl^£nYiJl  +  citAn^Tit^ ' 


j  *  1»  •••»  p  i  =  1 » 


,  s. 


(3.3.2) 


The  matrix  of  second  partial  derivatives  are  as  in  (3.1.6), 


(3.1.7),  and  (3.1.9)  with 


j  *  1 ,  . . . ,  p  i  =  1 ,  . . . ,  s.  (3.3.3) 

To  calculate  Ed^oj)),  Edj1^^))  and  E(l|2)(aj)) 
consider  any  of  the  s  stress  levels.  The  n^  observed  failures, 
Y^,  have  the  following  survival  function 

F<(y)  *  exp  -  (  ?  A..y  j  )  ,  y  >  0,  i  =  1,  ....  s.  (3.3.4) 
j=l  J 


Let  f^£  denote  the  number  of  failures  in  the  interval 
[xi  £_j,  ri£),  £  *  1,  ....  +  1  where  xiM  +1  =  ~.  Define 

Fi£  =  and  ^i£  =  ”  fri£*  and  Ui£k’  k  =  ^ . fi£ 

denote  the  failure  times  of  the  f^.  items  which  fail  in  this 

interval . 

Cohen  (1963)  shows  that 


E<fu>  * 


NFU  for  t  »  1 


£-1  c 


(N  -  l  —  )(FU  -  Fu.-,)  for  £  =  1 . M.  +  1 

k=1  Fik  (3.3.5) 


if  cM^  is  fixed  and  for  £  =  1,  ....  if  clM  is  random 


Now 


fi  a 

* EE  y,  wfu) 


E(fU)E(UttJk|Tu  <  Um  <  TU) 


E(fU). 


/Uu“jdy,yVl)  <*<>(  - 1  -  Fu.,) 


for  £  .=  1 ,  . . . ,  M.. 

where  E(f^  )  is  given  in  (3.3.5)  and  the  integral  must  be 
evaluated  numerically,  with  similar  expressions  for 


(3.3.6) 


E  (ff  UUktnUUk)  Mtl  E  (X  “ukl^W2)-  Usi"9  these 

k=l  K" 1 

expressions  when  a  fixed  number  of  items  are  removed  at  time 


fi£  a. 


these 


Tijy|  we  have 


E<T-i (“i) )  *  Ni  /  u  j  f  QtkX.  -u  k_1 )  exp  (  -  l  X-ku  k)  du 
1  J  1  Q  "k=l  *  1J  k=l 


i  «  a,  kfl 
Ci£  /  “  * 


L  akxiku 


ii  • 

exp(-  J  X1kU°k)du  ♦  j’ 


E(T|1)(“j))  *  Ni  /Q  u<,Jlnu(  j,  “kx1k‘,C'k‘1  J  • 

P  ak  ^i  a. 

exP(-  jiXiku  K)du  +  ^  cuTiJ£n(Tu) 


(3.3.7) 


‘  t-1  r4  (  "aj,l"u  (k-i  “kX'k“Clk’1  )exp  j1V°k 1 du 

ri£  Ti£  * 1  k  1 


(3.3.8) 


and 


E(Tj2)(aj)) 


Ni  ^  u  %nu)2  (^^akXiku  k_1)  exp  (-  Xiku  k)  du 


M, 


-  X  A*nu)2  (  j.  akxikuVl) 


*=1  ru  Tu 


k=l 

M. 


exp  (’  j,AfkuIlk)  du  +  j] 

i  —  1  *  • .  • ,  s  j  —  1 ,  . . « ,  p  • 
When  all  testing  stops  at  time  there  are 


(3.3.9) 


M.-1 


M, 


CM4  Ni  "  X  CU  "  X  fi£  items  removed  from  test.  Thus 
i  A“1  t“l 


H,-l  c 


EVr«.  c-i,  r1 

i  i  i-i 

And,  here, 

ai  B  Ok-1 

E(T1(aj))-ll1{/o  u  J  ckX1kU  k  ) 

exp  (  -  jt  xjku°k)  du  +  rfH)  t“^> 


TiMj  a . 


^i~^  Cit  xiMj  d  a.-l  p 

-  L  ~L  1  u  J  (  1  Viku  k  )  exp  ("  l  Aiku  k) du 

i«l  Fu  Ti£  k=l  *  1K  k=l  1K 


(3.3.10) 


Y  a,  ”1-’  F1M,  a. 

+  t  cutu  ■  i  cu  — -  , 

A-1  1X  ,x  £=]  u  jr  ™. 


U 


E(T<1)(aj))  . 


t1M,.  a 


in.  a.  p  .  n  a 

Ni{/0  U  J*nU  (Jj  akXiku  )  «pC  -  kf  S-ku  k)  du 

+  F1MiT1M1£r,TiM1}  '  L  ZT-  U  lMlu  j  Anu  (  £  Vik^  ^ 
11  1  *al  FU  TU  k=1 


«  exp  ( -  f 


a- 


M.-l 

l  a. 


and 


Jt  AikU  )  *  +  FiMi  TiM.  AnTiM.}  +  l=]  cUTij  *"V 

(3.3.11) 

E(T|2)(°j))  = 


Ni(/T,M,  u“j(lnu,J  l«p(-  j/ik-“k)  du 


Mj-1 
1  C, 


+  riM,  -  J,  F~Cri"i  ’IN/*"1!*/ 

*^1H  (X  D 

*  lH  '  “  J(tnu,!  ( j,  °kVV,)exp(-  ^  Xfku“k)du  J 


M4-l 

1  a. 


*  J,  eUTit(tnTU>’-  (3-3.12) 

Also  when  Is  random,  n^  Is  a  random  variable  with  mean 


(3.3.13) 


Mrl 

E(nj )  *  NFjjj  -  £  c,j£ 

1  1^  fa  u 


FiM, 


it 


Substituting  these  expectations  In  the  appropriate  places  In 
equations  (3.1.7)  to  (3.1.10)  yields  the  asymptotic  covariance 


matrix  of  (8j0! 


,  ctj)  for  this  censoring  scheme. 
3 


3.4  Estimation  of  Parameters  at  the  Usage  Stress 


A  A  A 

Let  ay  g,  ...»  be  the  maximum  likelihood  estimators 

j 


Of  Oj,  Bj0, 


>.  k ,,n,  •••»  3,-l  t  j  =  1,  •••»  P  obtained  from  an  accelerated 

life  test  as  described  in  section  3.1  to  3.3.  Let  k  denote  the 

a  a  a  J 

covariance  matrix  of  (Bj0.  ....  Bjk  ,  otj).  Recall  ^  is  of  the 

3 

form 


(3.4.1) 


^(J.ot) 

/v 

Is  the  vector  of  covariances  between  and  (&jq»  •••»  Bjk  ) 

(Bjo»  •••»  ejkJ  and  °oa  the  variance  of  °j-  Let  Ej  be  an 
estimator  of  Jj.  Let  Vy  denote  the  design  stress  of  the  system. 
For  sufficiently  large  sample  sizes  the  vector 

/V  A  A 

(Bjq.  ....  3jk  »  <Xj) 1 s  approximately  normal  with  mean  vector 
(8j0»  ....  aj)and  covariance  matrix  Ej.  Following 

J  ... 


where  ^  is  the  covariance  matrix  of  (Bj0. 


Thomas,  Bain  and  Antle  (1969)  we  recommend  sample  sizes  of  at 
least  a  hundred  at  each  stress  level. 

At  stress  Vu  the  maximum  likelihood  estimator  of  the  scale 
parameter  of  the  jth  components  time  to  failure  distribution, 
Aju»  *s 


\lu  ■  exp(io  Vj*(V)  *  J  =  1 . .  <3-4-2> 

For  sufficiently  large  sample  sizes  Xju  has  a  log  normal  distri¬ 
bution  with  mean  Xju  and  variance  aju  given  by 


°ju  *  <’•  ejl<wu>- ••••  ejk.<Vu»Ij;j<>.  6 j , ( Vu ) . 8jk  (V„))T. 

J  J 

(3.4.3) 

Hence  a  reduced  biased  estimator  of  is 

J  ^ 


\ju  =  \ju  exP{-°ju/2*’  j*1 .  P  (3.4.4) 

A  ^ 

where  cxju  is  obtained  by  replacing  Jjj  by  in  (3.4.3).  If 
L*.  were  known  the  mean  squared  error  of  X.  is  *L(exp(o?  )-l) 

J  J  J  U  J  u  J  U 

which  is  always  smaller  than  X?  (exp(2o?)-2  exp(o^,)+l ) ,  the 

a  J  U  J  U  J  u 

mean  squared  error  of  Xju.  An  asymptotic  (1  -  a)  x  100%  confi¬ 
dence  interval  for  X.  is  given  by 


(Xjuexp(-Zi_a/2°ju^  *  \juexp^Zl-a/2aju^  »  J 


*  1 


p.  (3.4.4) 


The  maximum  likelihood  estimator  of  the  j 


th 


cumulative  hazard  rate  at  stress  V  and  time  t,  A 

by 


components 
ju(t)  is  given 


V^'^V  . p-  <3-4-5> 


This  estimator  is  also  biased.  An  asymptotically  unbiased 
estimator  of  Aju(t)  is 

SJBU)  *  j  =  t  >  0 

(3.4.6) 

where 


oj(t)  =  (l,  e^oy,  ejk  (Vu),  tnt) 

?jd.  eji(y.  ••••  ejk  (V„).  *nt)T. 

J 


(3-4.7) 


A 

This  estimator  also  has  smaller  mean  squared  error  than  A-  (t). 

(t) 

Asymptotic  (1  -  a)  x  100%  confidence  Intervals  for  Aj^'  are 
given  by 

(Aju(t)exp(-Z1_a/2oj(t)),  Aju(t)exp(Z1_a/2a;j(t))t 

3  ■  1 . P  t  >  0.  (3.4.8) 

The  maximum  likelihood  estimator  of  the  components 
survival  function  at  time  t  and  stress  Vy  Is  given  by 

FjuU)  a  exP(“AjU(t))»  j  =  1 . (3.4.9) 

Approximate  (1  -  a)  x  100%  confidence  intervals  for  the  jth 
components  survival  function  at  time  t  and  stress  Vy  are  given 
by 

j  *  1,  ....  p  t  >  0  (3.4.10) 

Let  <  be  a  subset  of  1,  ...»  p  of  cardinality  k.  We  are 
Interested  in  obtaining  estimators  of 

ftht)  ■  irjt),  (3.4.H) 

the  survival  function  of  an  item  which  can  fail  only  from  the 
failure  of  components  indexed  by  elements  of  k.  When 
<=  {1,  ...,  p}  then  (3.4.1)  is  the  overall  system  survival 
function.  When  k  is  a  proper  subset  of  (1,  ....  p)  then  (3.4.1) 
represents  the  survival  function  of  a  system  which  has  been 
redesigned  so  that  those  components  indexed  by  <c  are  extremely 
reliable.  Clearly,  the  maximum  likelihood  estimator  of  F^(t) 
is 

^(t)  «  n  f.(t).  (3.4.12) 

JCK 

An  approximate  (1  -  a)  x  100%  confidence  interval  for  F^(t)  is 


(3.4.13) 


JCK 


J£K 


JU' 


1+M  al^k 

where  v  =  *■  - .  This  is  a  conservative  interval  in  the 


Y - 1 

following  sense.  From  (3.4.10) 


,  xl/k  exp(Z  oi(t))  a  exp(-Z  o^t)) 

(l-a)1/k=P(Fju(t)  rJ  <rju(t)<Fju(t)  rJ  )  for  je<. 

A  A 

Since  (FJU(t),  Fj,u(t))  are  asymptotically  independent  for  j^j', 

a  exp(l  oAt))_  a  exp(-Z  0j(t)) 

(l-a).p(Fju(t)  r  J<  F1u(t)<r1u(t)  r  J  .for  all  je«) 


JU'”7-  JU' 


exp(Zjj.j(t))  a  exp(-Z  o-.(t)) 

<p(nr  (t?  ^  ifjuWi" r  (t)  ^  )■ 

J£K  °  J  J£K  J 


3.5  Dependent  Risks 

In  sections  3.1  -  3.4  it  was  assumed  that  the  component 
lifetimes  were  independent.  This  assumption  may  be  relaxed  by 
considering  a  fatal  shock  model.  For  simplicity  we  shall 
illustrate  this  model  for  the  bivariate  case. 

Let  Uj ,  Ug,  be  Independent  Weibull  random  variables 
with  shape  parameters  a-j,  and  scale  parameters 

Aj(V,  BjK  A2(V,  j^),  and  A^(V,  fJ^)  ™  an  environment  charac 
terized  by  a  constant  application  of  a  stress  V.  Here  Uj 
represents  the  time  until  a  shock  destroys  the  first  component 
only,  U2  the  time  until  a  shock  destroys  the  second  component, 
and  U-|2  the  time  until  a  shock  destroys  both  components.  If 
(X'j ,  X2)  represent  the  component  lifetimes  then,  clearly, 

Xj  *  min(Uj,  U^)»  and  X2  =  min(U2»  U^)*  The  component  sur¬ 
vival  functions  are  not  Weibull,  but  are  given  by 


Fj(t;  V)  =  exp(-Aj(V,  B.)taj  -  A12(V,  612)t“12) • 

j  *  1,  2  t  >  0.  (3.5.1) 


An  accelerated  life  test  can  be  conducted  as  before  where 
now  the  parameters  of  interest  are  ,  02,  a^»  £] » 

The  results  of  section  3.4  can  be  used  to  obtain  estimators  of 
the  component  survival  function  under  normal  conditions. 

4.  EXAMPLE 

As  an  example  of  these  procedures  we  shall  consider  an 
example  given  in  Nelson  (1974a).  The  problem  is  to  analyze 
an  accelerated  life  test  conducted  on  Class-H  insulation 
systems  for  electric  motors.  There  are  three  possible  types 
of  insulation  failures  corresponding  to  distinct  parts  of  the 
insulation  system,  namely  Turn,  Phase,  and  Ground.  The  failure 
cause  is  determined  by  an  engineering  examination  of  the  failed 
motor. 

The  purpose  of  the  experiment  is  to  estimate  the  average 
life  of  such  insulation  systems  at  a  design  temperature  of 
180°  C.  A  median  life  of  20,000  hours  is  necessary  for  the 
satisfactory  performance  of  these  insulation  systems.  To  reduce 
test  time  and  cost  an  accelerated  life  test  was  conducted  at  4 
accelerated  temperatures,  namely,  190°  C,  220°  C,  2400  C,  and 
2600  C. 

The  accelerated  life  test  was  conducted  by  putting  10 
motors  on  test  at  each  of  the  4  stress  levels.  Motors  were  run 
until  they  failed,  then  the  cause  of  failure  was  found  and 
isolated  and  motors  were  run  until  a  second  failure  occurred. 

The  results  of  this  study  are  reported  in  Nelson  (1974a).  The 
data  followed  a  log^Q  normal  distribution  so  the  Weibull  theory 
results  do  not  apply. 

To  illustrate  the  results  of  the  previous  section  Nelson's 
example  is  reproduced  by  simulating  the  life  test  using  a 
Weibull  model  with  shape  parameter  1  for  each  failure  cause. 

The  shift  parameters  are  chosen  by  fitting  an  Arrhenius  Reaction 
Rate  model  to  the  estimated  component  medians  obtained  by  Nelson. 
The  model  is 


4760.6289  Turn  350.6594  Turn  238.5782  Phase  1315.0032  Ground 

1621.5481  Turn  2026.7441  Phase  1001.3643  Ground  90.0674  Turn 


(4.1) 


(V»  6.j)  =  6xp(BjQ  +  B j i © j i ( V ) )  s  j  =  1,  2,  3 

where  e^(V)  =  -1000/V  for  j  =  1,  2,  3  and  V  is  the  temperature 
in  degrees  absolute.  The  absolute  temperature  is  273.16  plus  the 
centigrade  temperature.  The  values  of  (Bj0*  Bj-j),  j  *  1,  2,  3 
are  as  follows: 

Table  4.2  True  Values  of  6^ 


Bo  Bi 

Turn  8.2607  8. 0106 
Phase  3.7748  6.1253 
Ground  13.0340  10.6487 


Twenty  Weibull  observations  were  generated  at  each  of  the  four 
stress  levels.  The  data  are  in  Table  4.1. 

A  Newton -Raphson  procedure  was  used  to  solve  the  likelihood 
equations.  The  integrals  in  (3.1.14),  (3.1.15),  (3.1.16),  and 
(3.1.17)  were  evaluated  using  a  repeated  seven  point  Gauss- 
leguerre  formula.  The  maximum  likelihood  estimates  are  as 
follows: 


TURN 

PHASE 

GROUND 


a  «  1.0099,  80  =  6.1363,  B]  =  6.9390 
a  =  .9993,  8q  =  3.5831,  =  6.0272 
a  =  1.0395,  8q  =  7.9788  B]  =  8.2679. 


The  estimated  covariance  matrices  are: 

/N 


^TURN  s 


^PHASE  B 


^GROUND  " 


/9.8252 

5.3254 

.1178 

(  5.3254 

3.1306 

.1283 

V  .1178 

.1283 

.0188 

/21.7658 

11.6281 

.2217 

11.6281 

6.6803 

.2442 

\  .2282 

.2442 

.0353 

A  8. 2378 

9.8173 

.1856 

1  9.8173 

5.6782 

.2033 

\  .1856 

.2033 

.0299 

) 


Using  these  estimates,  at  a  design  stress  of  180%  the 
estimates  of  the  probability  of  component  survival  at  a  mission 
time  of  20,000  hours  is  .1020  for  turn  failures,  .3023  for  phase 
failures,  and  .3574  for  ground  failures.  90%  confidence  intervals 
for  the  probability  of  component  survival  at  20,000  hours  and  a 
temperature  of  180°  C  are: 

TURN  -  (.0187,  .2701), 

PHASE  -  (.0749,  .5756), 

GROUND  -  (.1191,  .6080). 

Using  equations  (3.4.12)  and  (3.4.13)  the  maximum  likelihood 
estimate  and  90%  confidence  interval  for  system  reliability  at 
20,000  hours  and  a  temperature  of  1800  C  are  .0110  and  (.000027, 
.1402).  Similarly,  a  90%  confidence  interval  for  a  redesigned 
system  in  which  turn  failures  cannot  occur  is  (.0043,  .4024). 

We  note  that  the  above  confidence  intervals  are  suspect  due 
to  the  relatively  small  sample  sizes  and  are  provided  here  to 
only  to  illustrate  this  procedure. 
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